METHOD AND APPARATUS FOR NON-INVASIVE 
BLOOD CONSTITUENT MONITORING 

BACKGROUND 

The present invention is related to U.S. Patent Numbers 5,372,136 and 5,499,627 the text 
and drawings of which are incorporated herein by reference as if reproduced in full. 
THE FIELD OF INVENTION 

The present invention relates to improvements in the systems and methods for non- 
invasively measuring one or more biologic constituent concentration values. More particularly, 
the present invention relates to non-invasive spectrophotometric systems and methods for 
quantitatively and continuously monitoring the hematocrit and other blood parameters. 
THE PRIOR ART 

Modern medical practice utilizes a number of procedures and indicators to assess a 
patient's condition. One of these indicators is the patient's hematocrit. Hematocrit (often 
abbreviated as HCT) is the volume expressed as a percentage of the patient's blood which is 
occupied by red corpuscles, commonly referred to as red blood cells. The present invention is 
presented in the context of hematocrit. However, it is to be understood that the teachings of the 
present invention apply to any desired biologic constituent parameter. 

Medical professionals routinely desire to know the hematocrit of a patient. In order to 
determine hematocrit using any of the techniques available to date, it is necessary to draw a 
sample of blood by puncturing a vein or invading a capillary. Then, using widely accepted 
techniques, the sample of blood is subjected to either high-speed centriftige, cell counting, 
ultrasonic, conductometric or photometric methods of evaluating the sample of blood in a fixed 
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container. Prior U.S. Patent Number 5,372,136 indicates a system and methodology for 
determining the hematocrit non-invasively, without puncturing or invading the body, 
spectrophotometrically and continuously in a subject. The present invention relates to 
improvements upon the above cited system. 

Beyond the above referenced patent, others have suggested various means of noninvasive 
measurement of hematocrit. Specifically, Mendelson, U.S. Patent Number 5,277,181; Seeker, 
U.S. Patent 5,188,108; Gonatas, U.S. Patent Number 5,528,365; Ishikawa, U.S. Patent Number- 
5,522,388; Shiga, U.S. Patent Number 4,927,264; Tsuchiya, U.S. Patent Numbers 5,441,054, 
5,529,065, 5,517,987 and 5,477,051; and Chance, U.S. Patent Numbers 5,353,799, 5,402,778, 
and 5,673,701 have attempted to define means of directly measuring desired biologic constituents 
such as hematocrit. Even though the various patents indicate the need to utilize muhiple 
^wavelengths measured at different detection sites and/or the need to perform differential or 
ratiometric operations on the detected optical signal, all fail to isolate and resolve the individual 
and specific scattering and absorption coefficients of the desired constituent. At best they address 
only bulk attenuation coefficients and/or bulk diffusion constants of the scattering media while 
attempting to resolve such constraints as tissue nonhomogeneity. As an example, tissue may be 
considered to contain a bulk absorptive coefficient due to blood, collagen, water, fibers, bone, 
fingernail, etc. Hence, in order to determine the absorptive coefficient of the blood itself, the bulk 
value of the tissue per se must be prorated by the amounts of the above constituents. Secondly, 
the actual absorptive coefficient of the blood must then be decoupled or isolated fi*om its 
proration factor as well. 
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OBJECTS OF THE INVENTION 

Thus, it is an object of the present invention to provide an improvement in the systems and 
methods for the non-invasive (transcutaneous) and continuous determination of the blood . 
Hematocrit in living tissue. 

It is yet another object of the present invention to provide an improvement in the systems 
and methods for the non-invasive (transcutaneous) and continuous determination of the blood 
constituents, including glucose, bilirubin, cholesterol, tissue water, etc. in living tissue. 

It is another object of the present invention to provide a system and method and apparatus 
for the display of both immediate and/or continuous visual information regarding the HCT of the 
subject. 

It is yet another object of the present invention to provide a repeatable and 
reliable method and apparatus for the non-invasive determination of hematocrit transcutaneously 
and in real time even under varying physiological conditions. 

Still another object of the present invention is to provide a method and apparatus for the 
instantaneous determination of the bulk absorption coefficient of the scattering media. 

These and other objects and advantages of the invention will become more fiilly apparent 
from the description in the specification and claims, which follow. 
SUMMARY OF THE INVENTION 

In one aspect, the present invention accomplishes the transcutaneous, noninvasive, real- 
time and continuous measurement of the hematocrit and other blood constituents of the patient. 
That is, the electronic circuitry necessary is included to receive signals from a detector and to 
generate appropriate signals at various input sites as described in U.S. Patent Number 5,372,136. 
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Yet another aspect of the present invention is the ability to extract the blood absorption 
coefficient from the bulk tissue diffusion constant or the bulk absorption coefficient of the 
scattering media by requiring both physical and mathematical operations. 
BRIEF DESCRIPTION OF THE DRAWINGS 

Figures 1 and 1 A show a finger placed into a clam-shell type fixture constituting a 
receiving means for detector and emitter arrays operating in a transmission mode and the blood 
conduit which in the figures is the finger. 

Figures IB and IC are similar to Figure I A, but show the detector and emitter arrays 
operating in a reflectance mode. 

Figure ID is a schematic diagram for a mylar base with a detector, emitters and either a 
strain gage or a pressure transducer for inclusion in the clam-shell fixture. 

Figure IE is a schematic diagram for a detector emitter array using a single, moveable 
emitter in a transmission mode. 

Figure 2 shows actual patient data plot of ln(i) vs. d. 

Figure 3 illustrates actual patient data of the (di/dt)/i dependence on d. 

Figure 4 shows the blood absorption coefficient's dependence on hematocrit. 

Figure 5 indicates the nonlinear relationship between the V^/Vj-and pressure. 

Figure 6 shows the electrical circuit diagram of the piezo film / strain gage transducer 

means. 

Figure 7 is the plot of /(H) vs. measured Hematocrit. 

Figure 8 shows the instantaneous time derivatives of {^i/(?t)/i and d P/dt versus time 
during one cardiac pulse. 
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Figure 9 plots {di/dt)li versus d P/dt for a given human pulse at d^, d^ and d^ 
Figure 10 plots {di/dty{d P/dt) versus time during a single cardiac pulse cycle. 
Figure 1 1 is the circuit diagram of the pressure transducer means. 
Figure 12 plots versus at a fixed Hematocrit. 

Figure 13 gives the patient data of the new transcutaneous Hematocrit method and system 
plotted versus the measured Hematocrit standard. 
DETAILED DESCRIPTION OF THE INVENTION 

In a preferred embodiment of the invention, measurements are conducted using a modified 
version of the apparatus described in U.S. Patent Numbers 5,456,253 and 5,372,136, both of 
which are incorporated herein as if reproduced in fiill below. Both of these patents form part of 
the present disclosure. 

Thus, in a preferred embodiment, hematocrit is measured in living tissue located at some 
convenient location on the body, such as, an ear lobe, finger tip, nose or other accessible tissue 
sites. In a preferred embodiment the apparatus and signal manipulations described in U.S. Patent 
Number 5,372,136 are utilized to measure various optical parameters that will be described 
hereafter. The numbered components in Figures I, 1 A, IB, and IC are similar to the numbers in 
Figure 1 of U.S. Patent Number 5,456,253. 

In the present disclosure, Figure 1 shows the finger 7 of an individual placed into a clam- 
shell type fixture 6 wherein the optical and other physical measurements can be easily performed. 
The clam-shell type holder allows for adaptability to various finger sizes. However, other fixture 
methods such as Figures IB through IE, can be used to obtain similar physical data as using the 
clam-shell fixture. 



THEORETICAL BASIS OF THE SPECTROPHOTOMETRIC AND MATHEMATICAL 
ANALYSIS FOR TRANSCUTANEOUS HEMATOCRIT MEASUREMENT 

Non-invasive, transcutaneous hematocrit measurement using a spectroscopic method is 
described below: 
I. Introduction 

Earlier spectrophotometric techniques have fallen short of being able to fiilly characterize 
the individual blood absorbance coefficients. The following discussion demonstrates the method 
of decoupling, or isolating from the bulk tissue attenuation parameters (including the convoluted 
absorptive and scattering parameters) the individual blood absorptive constants. This unique 
method identifies, isolates and compartmentalizes each of the contributing biologic elements of 
the tissue media. This decoupling process can either isolate the blood absorbance of interest 
and/or eliminate the scattering contribution from the bulk media measurement. 

From photon diffusion analysis: 

SLjV(p)-c^'i!(p)='-m (1) 
dp' D 

where, 

1 

2(K + S) 
a = \/3K(K + S) 

K, = K iaJSAT+ (1 - SAT)a„) + (1 - H)K^ 
V 



(2) 

(3) 
(4) 

(5) 
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(6) 



V 



V) 



and where, 



a 



= Bulk attenuation coefficient of the tissue sample 



K 



= Bulk absorption coefficient of the tissue sample 



S 



= Bulk scattering coefficient of the tissue sample 



D = Diffiision constant 

Ki, = Macroscopic absorption coefficient for whole blood (WB) 

Sf, = Macroscopic transport - corrected scattering coefficient for WB 

K^. = Macroscopic absorption coefficient for plasma 

AT, = Macroscopic absorption coefficient for skin, & other non water/blood 
components 

= Macroscopic absorption coefficient for water 

V = Volume of a red blood cell (RBC) 

H = Hematocrit, volume fraction of RBCs to total blood volume 
SAT = Oxygen saturation % 

= Absorption cross - section of oxygenated RBCs 
o^, = Absorption cross - section of deoxygenated RBCs 
a, = Transport - corrected scattering cross-section of RBCs 
Xf, = Fractional volume of blood per total tissue volume 

= Fractional volume of skin, & non water/blood components per total tissue 



volume 
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= Fractional volume of water per total tissue volume 
^(p) = The photon density at a distance p 
^(p) ^ The source function, 
n. Analysis 

The light flux, or intensity, /, is given by / = Z?3y . When evaluated at p = d, one 

bp 

solution to equation (1) is: 

/ = A_s!i_ (8) 
e^"^- 1 

where A is a nontrivial function of the tissue scattering coefficient, S, the distance, d (if 

small), jnd the bulk attenuation coefficient, a. \iad» I, then (8) becomes: 

i^Ae"^ (9) 

where A _a 

[d" • (l-e-^«')] or (l/d^ + l/ad) for 0<n<2, 

where n is the power that d is raised to. 

Figure (2) shows the actual patient data plot of ln(i) vs. d, where a is determined 
directly from the slope of the line. 

The attenuation coefficient, a, is a bulk term which encompasses the attenuation 
measurement sensitivity to variations in skin color, presence of bone, callous, blood and 
water content, etc. In addition, a expresses the optical "path lengthening" effects of both 
the absorption and scattering characteristics of the tissue. Therefore, since a is a 
function of HCT and the intensity of the transmitted light can be measured, the HCT can 
be calculated by manipulation of the preceding relationships. 

Beginning with equation (9), the troublesome arid complex tissue function. A, 

8 



can be eliminated by taking the logarithm of (9) and differentiating with respect to the 
distance, d. Unfortunately the term A"^ is not known but changes with time as a result of a 
patient's cardiac cycle. Therefore, by differentiating with respect to time, this parameter 
becomes the time rate of change of blood volume which can be obtained through several 
methods described below. These time and distance derivatives may be performed in 
either order. 

[1] Taking the logarithm of (9) and differentiating with respect to the distance, rf, 

yields: 



a = a[ln ( / ) ] 
dd 



(10) 



Next the derivative of (10) with respect to time, t, gives: 



da 




(H) 



dt 



dt 



[2] Alternatively, first differentiate (9) with respect to time, t, to get: 



dt dx^ a dx, a dx^ at 



(12) 



When di dK, and 5i ^„ are negligible, and normalizing (12) by i yields: 

dx, a dx^ a 



di/dl = da (d - d„) , where d, -J, 2d 

Jlai J 



(13a) 



/ dl a f 

Figure 3 plainly demonstrates the oflfset term when the various graph lines are 
extrapolated to d = 0. The amount of ofiFset is shown along the y-axis. 

Next differentiate (13) with respect to distance, d, to eliminate that offset term to 



get: 



dd 



dt KdXj dt 



(14) 



Equations (3) - (7) are now used to extract the hematocrit from a. Squaring (3) and 
differentiating with respect to time results in: 



2a — = J 



dt^ ^ dl 



(15) 



Substituting the derivatives of (4) and (6) into (15) and rearranging: 



oa _ J 
a7 ~ 2a 



dl 



dl 



dl 



\ dl dt . 



(16) 



At 805 nm 



« S so that (16) can be simplified to: 
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dt 2a dt 



(17) 



By using the 805 nm wavelength the red blood cell absorption cross-section constants 
are equal, = a„, and /T^ is negligible. The hematocrit can then be determined directly 



from Kj, as (5) simplifies to: 



(17a) 



Figure 4 shows the linearity ofK^fH). 

IfK^ » KS^, where S is approximately 1.0 / mm in human tissue, then solving (17) 
for Kb and substituting into ( 1 7a) gives: 
da 



To rewrite in terms of measurable intensity, i, (10) and (14) are substituted into (18) 
to obtain: 



IfK^ is not » KSt. then substituting (5) and (7) into (17a) and rearranging terms 



W a[ln(0] d [{diidt)/i] 
3a dd dd 

a ■ 



(19) 




yields: 




(I8a) 
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Alternatively from (13a) : H » a • (dildtli>i 

(d-d J • 



(18b) 



Equation (18a) indicates a small nonlinearity in H may occur based on the magnitude 
of K for a given individual. 

It should be reiterated that the change in received intensity with time is a result of the 
change in normalized blood volume resulting from the cardiac cycle itself as blood pulses 
through the examined tissue. As the intensity of the received light is measured, its time rate 
of change can be calculated. The change with distance can be determined by placing multiple 
emitters (such as 1-4 in Figure 1 A) and/or multiple detectors such that multiple thicknesses 
of tissue and hence, lengths of tissue are penetrated. 

To examine 32^ further, the following can be defined for the illuminated tissue: 



Vf^ = Volume of blood, 



Volume of water, and 



= Volume of skin, tissue and other non-water or blood components. 



By definition. 



V, 



(20) 



differentiating (20) with respect to time gives: 



dt 



(21) 



\ 



Since dl^ « dYk and « + V„ (21) simplifies to: 
dt dt 
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dt (22) 

It is emphasized that a is a function of the bulk absorption and scattering coefficients, K and 
S, as well as hematocrit, H. Further, that K and S are functions of the fi-actional volumes of 
each constituent, X,_ and X„. which must be used to prorate the individual absorption and 
scattering coefficients, K^. K„ and S.. Therefore, the transducer system must be 
responsive not only to a change in volume (A^O due to the influx of the blood, but must also 
be responsive to the normalized change in volume of blood, normalized to the total volume of 
the finger (f^or tissue being measured. 

For Reflectance (R) measurements in homogenous tissue: 

R = Ae"*", where A « (1/H + 1/ar), where r is the radial distance, 

and 

iiE = a' (r/(a^ + ar)) 
R 

However, for tissue, which is typically non-homogeneous with a dermal and subcutaneous 
layer, the reflectance will not be a trivial function but can be described as approximately: 

R=[(C, + C2)exp(-C3-r)]/i' 
Where C, and C, are inter-related photon flux densities between the dermal layer 12 and the 
subcutaneous layer, 12a (see Figures IC and IE). Likewise, C, is a strong function of z„ Zj, 
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cCi, and i.e., the thickness of the dermis or dermal layer 12, subcutaneous layer 12a, and 
their respective a's. 

Since Cj* is a function of the inter-related photon flux densities Cj and/or C2 and if 

Xb\ does not equal Xb 2, then the slope C3' will not be nulled out by the Xb' monitors 

mentioned. Therefore, Xbj must be greater than Xb/. Then the pressure or piezo monitors 

will compensate correctly. The circular pressure balloon is ideal for not only sensing the 

change in a pressure, but also providing. a pressure against the dermis causing Xb/ to be 

small. However, recognizing that the penetration depth of the 800 nm light typically extends 

through dermal layer 12 into the deep tissue, subcutaneous layer 12a, a different wavelength 

selection is appropriate. Thusly, when the photons only penetrate into the dermal layer 12, 

C3' will only be a function of z, and Those selected wavelengths, as mentioned in U.S. 

Patent No. 5,372,136, would be the green (570-595 nm) wavelength and 1300 nm 

wavelength. The green wavelengths are used as the hematocrit bearing wavelength and the 

1300 nm wavelength is used as the non-hematocrit bearing, or reference wavelength. That is, 

for reflectance measurements the green (Gr)-1300 wavelength pair would give the hematocrit 

information as: 

AGr/Gr • _a^_ =/(HCT) 
A1300/1300 ai3oo 

in. Methods of ^ measurement 

dt 

dXJdi can be measured and compensated for through the use of a number of different 
methods - (a) a pressure transducer, (b) a strain transducer such as piezo electric film or 
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strain gage, (e) a difiFerent wavelength of light, such as 1300 nm, which also holds d^^dt 
information, but holds little hematocrit information, or (d) other transducers. The individual 
methods of obtaining hX^t are addressed below. 
A Pressure Transducer Measurement of ^ 

a/ 

Consider a pressure transducer system 36 with a gas filled bladder 38 surrounding a 
finger tip 10 of a patient contained within a fixed volume clam shell fixture 6, see Figures 1, 
1 A- ID. The same derivations, equations, and results would apply to any other body 
appendage or tissue that could be contacted such that a change in the tissue volume would 
change the pressure of the contacted pressure transducer system. For a finger note: 

y ^ = J/ + K, (23) 

where 

yciom - Clam-shell fixture volume 
V^^ - Bladder system volume 
Vf = Finger volume 
Also LVf= -t^V^ The system will have a bulk modulus of elasticity, P, such that: 

Substituting (23) into (24) results in: 



V 



da"* 1 I :— 
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Since A K- = then from (25) we have: 

• dJC^jKi^^^ (25a) 

As stated above, p is a constant of the pressure transducer system. However, an 
empirical solution for HjIsz. _ i 



was 



found to have a nonlinear relation to the pressure of 



the transducer system. For a given clam shell - pressure transducer embodiment a 



polynomial, F(p), can accurately describe 



dam ^ ^ 



, see Figure 5 



B. Strain Transducer (Strain Gage/Piczo Electric Film) Measurement of ^ 

. dt 

Again it is assumed that Af^^ = AP>. and that the finger changes volume only by a change 
in diameter. A strain gage or piezo electric film is secured tightly around the finger (again any 
applicable body appendage or tissue would apply) such that a change in diameter would produce 
strain in the transducer. Specifically assuming a cylindrical finger: 



dt dt dt dt 

Normalizing with respect to Vj, yields: 

a V^, ^zr rdt (27) 



A change in the length of the transducer element is related to a change in finger 
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radius by AZ, = 2KAr, therefore: 

m,-2Mm=2rO) (28) 

dt L, 

where y(t) = dL/dt is the rate of change in the strain as a function of time. For a strain 
L 

gage this value can be measured from an appropriate electrical circuit, see Figure 6, as it is 
proportional to the rate of change in the gage resistance. 

For a piezo electric film the voltage produced is proportional to the strain, therefore: 
2^ = _2_ dvm (29) 

dt g„v dt 

where, gj, is the piezoelectric coefficient for the stretch axis, t is the film thickness and v(t) is the 

open-circuit output voltage. 

C. 1300 nm Light Measurement of 

dt 

The selection of the 1300 nm wavelength is based on criteria established in U.S. 
Patent Number 5,372, 136. The approach here is not to solve for dXJd^ and substitute into (19) 
but to ratiometrically eliminate dX^t. In the case of the 1300 nm reference wavelength, the 
assumptions following equation (12) are no longer valid; i.e., dX/dt and dXJdt are not negligible, 
since water absorption at 1300 nm is so large. Hence, for the 1300 nm equations (13), (14) and 













3 


\dt) 


11 


~2a 



{(2 K^S)K,^KS,fjf-^-l2K^S)K,^ KS, ^ + 



(30) 
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where,, a, and the bulk and material specific K. and S are wavelength (A) dependent. 
Recalling that, Xi, + X, + X^= \,hy definition, and that: 

dt dt dt (31) 
By substituting (3 1) into (30) and noting that K^,j = ^^j. the following is obtained: 



f— 1 




3 




\i 


~ 2a 



(32) 



Since, » (32) becomes: 
dt dt 



(da 
dt. 



ij 



(33) 



Therefore, to eliminate dXh and solve for the hematocrit, (17) is divided by (33) 

dt 



yielding: 



{daldt\ _ao.Sii- 



(34) 



Since .S, and are well behaved and known (let K,j /S^ = G) in human tissue and 
the ratio Kis. is a function ofH, then rearranging (34) gives: 

■ ' N 



(35) 



dt. 
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Where da can be measured using (1 1) or (14). See Figure 7 for y(H). 
di 

D. Other dX^/dt measurements such as doppler, ultrasonic, electrical conductivity, magnetic 
permeability and other techniques have similar derivations. The important consideration is that 
dXy^t is a normalized time varying quantity. 
IV. Analytical implementation 

If hematocrit is constant over a given time interval, averaging can eliminate 
system noise whose frequency components have corresponding periods much shorter than the 
interval. In addition, by observing the data variance during the interval it may be determined that 
the data is invalid. In the present system, the data acquisition rate is approximately 1000 data 
samples per second. This means that within a typical human pulse about 1000 samples of data are 
available for appropriate numerical analysis, averaging and qualification. Recognizing that both 
the intensity of light and the pressure in the transducer system are changing in time during the 
influx of blood is of great importance. Since the parametric relationship ofda/dl as a function of 
^?/dt (where P is pressure) during the cardiac cycle should be linear, a multiplicity of data points 
facilitate qualification of the signal for accuracy and linearity. Whereas, prior techniques involving 
only the peak and valley values of the cardiac cycle require numerous pulses to qualify the data 
set. See Figures 8, 9 and 10. 

Fig. 8 shows di/dt/i as well as dP/dt verses time during the cardiac pulse - it is a pulse 
showing « 200+ data samples during the pulse. 

Fig. 9 shows (di/dt)Ai vs dP/dt showing that within one cardiac pulse 200 plus data samples 
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are linearly related, i.e. trace up out of the "0" origin up to a maximum value and then back down 
toward the origin again. 

Fig. 10 shows da/dt/dP/dt versus time during one single cardiac pulse with 200 plus 
samples of data from time 15-45 giving a value of about 4.5 thousandths. The data can then be 
averaged, as if 200+ individuals pulse (max-min) values were actually taken as present day 
oxymeters do. 

A. Homogeneity 

Since the above derivations are based on the assumption of tissue homogeneity 
(Le.,dXf,/dt = dXg,2/d t , A^=^ Aj, aA/aXb = aA^/aXb, a, = a^, etc.), high-speed, single-pulse, 
multiple parameter sampling allows for mathematical qualification of homogeneity, by requiring 
linearity of ln(i) vs. d and {dildt)li vs. d Under these constraints and when qualified as 
homogeneous, (da/dt)/(dP/dt) also may be assumed to be linear over the entire pulse contour. 
Finally, both a and da/dt must also be linear, fiirther assuring homogeneity in X^^ and in dXy/dt 

B. Circuitry 

See U.S. Patent Number 5,372,136 for the operational circuitry description, which 
allows for high speed sampling of the optical intensities. See Figures 6 and 10 for similar circuitry 
considerations for sampling of pressure, peizo, and strain-gage measurements. 

The circuitry shown and discussed in US Patent Number 5,372, 136 is programmable by 
conventional techniques to solve and implement the equations and calculations presented in this 
application. Figure 6 shows a piezo transducer circuit having a transducer 50 connected to a 
series of operational amplifiers, resistors and capacitors in accordance with the figure. The circuit 
terminates in an analog output 52 for connection to the "E" connection shown in the middle left 
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side of Fig. 9D in U.S. Patent Number 5,372,136. Figure 1 1, on the other hand, shows a pressure 
transducer circuit having a pressure transducer made 62 connected to a series of operational 
amplifiers, a capacitor, resistors and variable resistors as shown in the figure. The circuit 
terminates in an analog output also connected to the aforementioned "E" connection. 

Referring more specifically to FIG. 6, a crystal oscillator is connected to ground and to the 
non-inverting input of a first operational amplifier, which may be an LM158. The non-inverting 
input of the first operational amplifier is connected to ground by a .047 |iF capacitor C3. The first 
operational amplifier's feedback path to its inverting input includes a 470 K resistor R8. The first 
operational amplifier is suitably biased at the junction of a 220 Q resistor R7 and a 150 ^F 
capacitor C4 that are connected between VCC and ground. 

A second operational amplifier, which may also be an LM 158, receives the output of the 
first operational amplifier at its inverting input via a 10 KQ resistor R5. The second operational 
amplifier's non-inverting input is connected to several locations: 

to a voltage VB51, which may be 4.096 volts, through a 10 KD resistor R2; 

to a middle node of a voltage divider, the voltages divider extending between the 

non-inverting input of the first operational amplifier via a 10 Mfl resistor R4 to the middle 

node, and via a 10 KQ resistor Rl to ground; 

to the inverting input of the first operational amplifier via a 10 KQ resistor R9; and 

to ground via a 220 |iF capacitor C5. 
The second operational amplifier's feedback path to its inverting input includes a parallel 
arrangement of a 0. 1 )iF capacitor C2 and a 47 KQ resistor R6. The second operational amplifier 
drives the A/D output 52 via a 10 KQ resistor R3, the output connected to ground via a 1 ^F 
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capacitor CI. 

Of course, the particular choice, arrangement and values of components shown in FIG. 6 
may be varied while still remaining within the scope of the invention. 

Referring now to FIG. 10, first through fourth operational amplifiers, which may be 
LM348s, are illustrated. The operational amplifiers are powered and biased by voltages VCC and 
VEE. 

The first operational amplifier's non-inverting input is set to a value determined by the tap 
setting of a 1 KQ adjustable resistor R2 that extends between VCC and VEE. The DAC input 
drives the first operational amplifier's inverting input via a 1 KQ resistor Rl. The first operational 
amplifier's feedback path includes a 50 KQ adjustable resistor R4. The first operational amplifier 
drives the second operational amplifier's inverting input through an 1 1 KQ resistor R3. The 
feedback path to the inverting input of the second operational amplifier includes a 100 Q resistor 
R5. 

A transducer 62, which may include a Motorola MPX20100P, has opposite terminals that 
drive the non-inverting inputs of the second and third operational amplifiers, respectively. The 
other two opposite terminals of the transducer are connected to VCC and ground, respectively. 

The second operational amplifier drives the inverting input of the third operational 
amplifier via a 750 Q resistor R6. The third operational amplifier's feedback path to its inverting 
input includes a parallel arrangement of a 93.1 KQ resistor RIO and a .001 jiF capacitor CI. 

The third operational amplifier drives the non-inverting input of the fourth operational 
amplifier via a 1 KQ resistor R7. The inverting input of the fourth operational amplifier is 
connected to ground via a 1 KQ resistor R8. The feedback path to the inverting input of the 
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fourth operational amplifier includes a 50 KQ adjustable resistor R9. The fourth operational 
amplifier drives the output of the FIG. 1 1 circuit. 

Of course, the particular choice, arrangement and values of components shown in FIG. 10 
may be varied while still remaining within the scope of the invention. 
C. Preferred Embodiment 

Physical embodiments as shown in Figure I include the optical array, pressure 
transducer/balloon system and clam-shell fixture. Requisites of the preferred embodiment include 
a holder for the finger (or other tissue) such as seen in Figures 1 and 1 A and IB. This clam-shell 
fixture not only secures the tissue but also the optical array, and transducer system. 

Figure ID is a schematic diagram for a mylar base member 38 that is shaped generally like 
a cross. As oriented in Figure ID, vertically extending portion 52 crosses with a horizontally 
extending portion 54 to yield top leg 56, bottom leg 58, and side legs 60, 62. In use, a finger 7 lies 
along the longitudinally extending portion 52 with the finger tip placed on the top leg 56 to 
properly cover the arrangement of LED's 32 and photodetector 34, which are arranged like those 
on Figures lA - IC. A piezoelectric pressure transducer or strain gage 66 spans the horizontally 
extending portion 54 fi-om near the tip of side leg 60 to the tip of side leg 62. In this orientation, 
the transducer or gage may be wrapped around the finger 7 for use in measurements. 

The optical array 30, seen in Figure ID, shows the arrangement of multiple LED's 32 
spaced at known separation distances from the detector 34. This array provides for the 
instantaneous distance, or "tf derivative, by the transmission mode shown in Figure 1 A or in 
reflectance modes shown in Figures IB and IC. However, as shown in Figure IE, a single LED 
42 swept across the finger 7 or tissue surface 9 with a stepper motor 44 would provide a d 
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derivative as would a cantilevered clam-shell with an angular measurement device. In any case, d 
must be known and/or fixed. Also the detectors and emitters may be placed anywhere about the 
finger. 

The pressure/balloon, strain gage, or peizo transducer system incorporated within the 
clam-shell fixture ( see Section III, A, B, C and Figure 1 A) provides the contact surface area 
needed to define the dX^/dt 

High-speed sampling provides for a closer approximation of the instantaneous time, t, 
derivative, didi, as opposed to peak- valley values, see Figure 8. Therefore, the above 
embodiments allow for the direct measurement of In (i) at ^4 ^2. and d4 cotemporaneously, 
thereby determining the actual a of the sampled tissue. Likewise (di/dt)/i can be directly 
measured at dj d2„ dj and d^ cotemporaneously during the pulse which determines the 
instantaneous 3a/3t. 

The above mentioned optical array can be utilized transmissively and/or reflectively 
provided the separation distance between the detector and first emitter (d^) is greater than 3mm. 
D. Choice of Non-Ionizing Wavelengths 

Since hematocrit is an example of the desired biological constituent concentration value of 
interest, selection criteria of the preferred wavelength must include an understanding of equation 
(5). That is, a wavelength whose coefficients are small compared to and which are 

also insensitive to oxygen saturation status must be selected. Such wavelengths include 805 nm, 
590 nm, 569 nm and other isobestic wavelengths with negligible water absorption. While non- 
isobestic wavelengths, with small water absorption, could fianction, a second wavelength is 
needed to null out the oxygen saturation effects. 
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If the desired biologic constituent value of interest is the blood glucose, bilirubin, 
cholesterol or other parameters, then a second wavelength must be chosen. The first wavelength, 
805 nm, is used to measure the hematocrit, H, after which a Kp^os (the absorbance of plasma at X 
= 805 nm) can be determined. Then, knowing the H, a second wavelength, 570 nm, is chosen 
where Kp57o is less than Kpgoj, Similarly, if the first wavelength used to measure the H and the 
reference glucose, (glucose) is 570 nm, the second wavelength, 1060 nm, is chosen where 
Kp57o is much less than Kpio^o. In the case of bilirubin, the first wavelength used to measure the H 
and the reference bilirubin, Kp (bilirubin), is 570 nm, the second wavelength, 440 nm, is then 
chosen when Kp57o is much less than Y^^. The selection of these above mentioned wavelengths 
therefore assures uniqueness for the measurement of the desired biologic constituent. 

Additionally for glucose determination, recall that the 1300 nm wavelength is not 
hematocrit or hemoglobin dependent but will be glucose sensitive. This is primarily due to the 
dependence of the scattering coefficient on the difference between the index of refraction of pure 
water and glucose, i.e.: recall 

iSbg = H (I'H) a,8 (from equation 7) where: 

o.« = 87iV(Tl*8-l)'-bVX^ 
where r]\ = index of refraction of the RBC hemoglobin at 800 nm relative to plasma (the 
plasma index of refraction), and, 

5^,13 ^ (1-^ (^so from equation 7) where: 

o.i3 = 87i^V(Tl\,-l)'-bVA.^ 
r|'i3 = the index of refraction of glucose at 1300 nm relative to r|o. 
Therefore, the 8 13 ratio has both hematocrit and glucose information. Whereas the ttg-aVAP 
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(equation 18a) ratio has only hematocrit information. Therefore the differential combination of 
those ratios will be a strong fUnction of glucose only. 
E. Improved Accuracy Pulse Oximeter Device 

The accuracy of present day pulse oximeters suffers from 4 major problems: 
tissue perfusion (low and low dXy^t), d dependence (varying finger sizes), tissue 
nonhomogeneity ( the tissue penetration depth for 660 nm light is not the same as for 940 
nm light), and dependence (see equation (5)). 

All of the above mentioned deficiencies in pulse oximetry can be eliminated by 

understanding equation (13). Equation (13) indicates an "offset term", 1 3A . 

A dX, 

Hence, while merely dividing (A///);^j by (A/Z/^j mitigates the effect ofdXJdi, the ds do 
not completely cancel, thereby yielding the above mentioned problems. To improve pulse 
oximeter accuracy, a derivative is needed as in (14), which eliminates the "oflFset term". Hence, 
the ratio of {daldt\J{da/it)^ results in no H, d, or dependence and the use of the multiple 
LED array and high-speed sampling as mentioned in section IV qualifies the tissue as 
homogeneous. 

V. A Simplified Two Step Approach 
(A) Determination of H and X/, 

The bulk attenuation coefficient, a, can be easily measured with the optical array, at 805 
nm, utilizing equation (10) and as described in Section IV(C). Notice that at 805 nm, a is a 
strong function of H and Xf, since K^g K^^ are small, see Figure 12. 

Therefore, by knowing itself, H can be determined. X^, itself can be determined using a 
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strain gage in the following two step approach. Step One, measure the strain gage resistance 
when the finger is made bloodless, by squeezing finger, such as with a stepper motor. Step Two, 
measure the strain gage resistance when the finger is blood filled, for example by suction. . 
Mathematically, at 805 nm and when K„ are small, equation (3) is approximated by: 

« 3KS (36) 

or 

a^=^3[/:^J[5^,-f5^J (37) 
Substituting (5) and (7) into (37) yields: 




WithJifi and a measured and known, and with the o*s and S,. X, approximately constant, H 
can be solved with a quadratic formula or a polynomial fit. 
The strain gage determination of Xf, is as follows: 

Let = the volume of a bloodless finger. Let = the volume of blood filled 
finger, and again considering the finger as a cylinder: 

V,^ 7rPz=K,+ (39) 
i(B}z^V,^V,^V^ (40) 

and 




(41) 
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From equation (20) 




(42) 



Substituting (41) into (42): 




(43) 



Where the strain gage resistances are proportional to the radius, r and R, of the 
finger. 

(B) Determination of tissu e- water content X... 

Choosing the wavelength of 1300 nm. where K, and K„ are significant, the tissue 
water content, ;sr,„.can be determined. Recall that l-X^-X^ = X, and substituting into 
(3) yields: 



With a, J.. and determined and because iiTi, K,. K^, andS, are known coefficient values 

at 1300 nm, is solved with either a quadratic formula or a polynomial fit. 

RESULTS 



apparatus and the application of Equation 19 on numerous patients with a correlation of 



= 3({/:, - K,}X, + {K^ - K,}X^ + KJ 



(44) 



Figure 13 demonstrates preliminary results with 30 patients the application of the method and 



0.96. 
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As implied throughout, those skilled in the art will also appreciate that the methods for 
determining blood hematocrit values within the scope of the present invention may be adapted 
for determining other non-hematocrit biologic constituent values such as glucose, bilirubin, 
cholesterol, tissue water, etc. 

The present invention may be embodied in other specific forms without departing from its 
spirit or essential characteristics. While the foregoing described embodiments are to be 
considered in all respects only as illustrative of the claimed invention, they are not intended to 
restrict the scope of the claims. The scope of the invention is, therefore, indicated by the 
following appended claims rather than by the foregoing description. All changes within the 
meaning and range of equivalency of the claims are to be embraced within their scope. 
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